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We compute graviton scattering amplitudes in M-theory using Feynman rules for a 
scalar particle coupled to gravity in eleven dimensions. The processes that we consider 
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^ describe the single graviton exchange and the double graviton exchange, that in M(atrix) 

t^. ' theory correspond to the v 4 /r r - and t> 6 /r 14 -term respectively. We further show that the 
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v 6 /r 14 -term appearing in M(atrix) theory at two-loops can be obtained from the covariant 

eleven-dimensional four-graviton amplitude. Finally, we calculate the w 8 /r 18 -term appear- 
ing at two-loops in M(atrix) theory. It has been previously conjectured that this term is 



related to a four graviton scattering amplitude involving the i? 4 -vertex of M-theory. 
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1. Introduction 

M-theory is the 'quantum version' of eleven-dimensional supergravity. The low energy 
degrees of freedom and interactions of this theory are those of general relativity However, 
if one wants to calculate e.g. radiative corrections to graviton scattering amplitudes one 
has to know how to describe the short distance degrees of freedom of M-theory. M(atrix) 
theory [|1] is a proposal that in the infinite momentum frame these degrees of freedom can 
be described in terms of Dirichlet zero branes. Graviton scattering amplitudes that are 
usually hard to compute (or even infinite) using standard quantum field theory techniques 
can be calculated in terms of a simple quantum mechanical model. So for example, the 
f 4 /r 7 -term appearing at one- loop in M(atrix) theory has an interpretation as a single 
graviton exchange diagram in eleven dimensions, while the w 6 /r 14 -term appearing at two- 
loops in M(atrix) theory describes a correction coming from general relativity 0. An 
example of a quantum gravity correction is the i> 8 /r 18 -term appearing at two-loops in 
M(atrix) theory 0. 

In this note we would like to consider these three terms in the effective potential 
between two gravitons in some detail. In section 2 we will compute the v 4 /r 7 - and w 6 /r 14 - 
terms using Feynman rules for supergravity. In section 3 we will show that the relativistic 
correction can be calculated using the covariant four-graviton amplitude in eleven dimen- 
sions ||. In section 4 we will make some remarks about quantum gravity corrections. 
The appendix contains more explicitly some of our calculations. 

2. Single Graviton Exchange 

In this section we would like to consider the scattering of two gravitons in eleven 
dimensions. The Feynman rules that we need to do this calculation were derived a long 
time ago by Feynman |J and De Witt J7| (see also ||). The vertices involving gravitons 
are rather complicated so that we will address the problem using a simpler model that 
captures the important properties of the calculation. First, we will be assuming that one 
of the gravitons involved in the process is heavy so that it can be treated as the source 
for the gravitational field in which the other graviton is scattered. Therefore, we will be 
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+ u hhifip" , (2.3) 



considering a problem of potential scattering. Our second assumption is that we can treat 
the 'probe graviton' as a scalar particle. The vertices describing the coupling of a scalar 
to gravity are simpler than vertices involving gravitons (see 0). 

We consider the eleven-dimensionalu Einstein action coupled to a massless scalar field 

S = J d n xV=g (JrR + Ig^tpd"?) , (2.1) 

and we expand the metric around a flat Minkowski background 

9iiv=Tivw + hp V . (2.2) 

Here h^ v is a small fluctuation that describes the source graviton field. The action for the 
'probe graviton' ip then becomes 

1 1 

2^+2' 

where we have indicated explicitly the cubic vertex and only symbolically the quartic vertex 
whose explicit form can be found in 0. The source graviton has vanishing transverse 
velocity and it moves with the speed of light in the eleven direction [|| . It is described by 
the Aichelburg-Sexl metric whose only non- vanishing component is /i__ 

h - 15N « X (24) 

where p s _ = N s /R is the momentum of the source graviton and M p i is the Planck mass 
that is related to the gravitational constant by k^ = 167r 5 /M^. The action takes then 
the form 

S v = ^ J d xl x (-<pn<p + h ++ d + (pd + ip + ...). (2.5) 

In the following we will consider the single graviton exchange process in the potential 
scattering description. The light-like compactification which is described in terms of a 
non-relativistic Schrodinger equation is relevant for finite iV M(atrix) theory ||10|1 . The 



space-like compactification is described in terms of a relativistic Klein-Gordon equation 
which leads to additional terms in the effective potential that are of no relevance for finite 
N. We will further generalize these results to the double graviton exchange describing the 
relativistic corrections. 



We are using the signature ( — ,+,...,+). 



2.1. One Graviton Exchange in the Light-like Compactification 

In the light-like compactification the equation of motion for (p that follows from (|2? 



is 



I d 2 d \ I 

+ i \(p = V<p with V = - — /i__<9 r 2 . (2.6) 



2m dx 2 dr J ' 2m 

Here m = N p /R, r = x + /2 and we have used that f is in a state of definite p-. Equation 
( |2.6|) has the form of a non-relativistic Schrodinger equation for a particle with mass m in 
an external potential V. 

The S-matrix can be obtained using the standard quantum mechanics techniques. The 
free propagator for this massive particle is 



K (x f , Tf, X h Ti ) = 0(Tf - n) J _0 e ^"/-- )-^(r,-r, 



(2.7) 



where (x, p) denotes the nine-dimensional transverse space and E = p 1 /2m is the light- 
cone energy. The plane wave solutions to the equation of motion have a non-relativistic 
normalization 

To first order the S-matrix is given by 

5 <» = -i J fU^TfW*,, Tf, 3, T)W, T)K B (Z, r; ft, n) W „(x„ n) d ,,^ id r. 

(2.9) 
After evaluating the different integrals we geta 

(1) _ _N £ N IL _v^_ J_ 

where k = pj — Pi is the nine-dimensional momentum transfer. The effective potential is 
obtained by Fourier-tranforming 

(1) __i5Ayv^ 

V 'ff~ 16 i? 3 M p V r ' { } 



This is the result that was computed in p| , [ \ I. ] and 



Here and in the following we omit the factor 2m5{Ef — Ei) that expresses energy conservation 
from all the scattering amplitudes. 



2.2. One Graviton Exchange in the Space-like Compactification 

The previous calculation can be repeated for a space-like compactification. In this 
case the equation of motion for ip has the form of a relativistic Klein-Gordon equation: 

d 2 d 2 



dx 2 dx\ 



+ m z )ip = V(p with V = -h—(d X0 -m) z , (2.12) 



where m is the eleven-component of the momentum. The solutions to the free equation of 
motion are now plane- waves with a relativistic normalization 

**> = <wr*im eW "> (2 ' 13) 

where E 2 = p 1 + m 2 . 

The S-matrix to first order is now 

S (1) = -i f \p* oat (y)V(y)ipin(y)d 10 y. (2.14) 

After performing the integration and a Fourier transformation we obtain 



r(1) _ Nph—iVT^-i) 2 



V}ti = --£ - vv ; (2.15) 

eff R 2 VT^J 2 

This result agrees with . While in the light-like compactification the S-matrix contained 
only a i> 4 -term Q2.ll ), here we obtain a series in terms of v. The difference is due to the 



fact that in the second calculation (p obeys relativistic kinematics, while in the light-like 
compactification tp obeys a non-relativistic Schodinger equation. These additional terms 
have no relevance for finite iV M(atrix) theory and are subleading in the large iV expansion. 

3. General Relativity Corrections 

3.1. Double Graviton Exchange in the Light-like Compactification 

Corrections coming from general relativity are described by the diagonal term in 
the table appearing in 0. They have integer powers in k%i. The simplest example of a 
correction of this type is the i> 6 /r 14 -term appearing in M(atrix) theory at two-loops. It can 
be obtained from the double graviton exchange or equivalently a second order scattering 



process. Here we will evaluate the contribution involving two cubic vertices of the type 
( ^6|) . To second order the S-matrix described by ( |2.6| ) is given by 



S {2) = - ¥>*out(xfi T f) K o(xf, Tf\ x', t')V(x', t')K (x', r'\ x,r) 

J (3-1) 

V(x, t)Kq(x, t; x*i, Ti)(pi n (xi, Ti)dxidxdrdx dr' dx f . 

We can evaluate the x*i, Xf, x and x' integrals exactly. From the r, r' integrals we get the 
contribution 

/+oo p+oo _ 2 

dre ir(E f -^) / dte^ E f~^\ (3.2) 

-oo JO 

where we have transformed to a new integration variable t = t—t'. The first integral again 
gives energy conservation but the second integral does not converge. We can evaluate it 
using a trick of Feynman and Hibbs [l!|. Expression (3J2) can be written as 



2^(E,-g,)Bm g/ _ E + fe . (3.3) 

where E = p 2 /2m. The S-matrix to second order takes then the form 

5 (3) = ^^! l im /^^ I 1 (3 4) 

WMlf 32 e-0 J (2tt) 10 \p f - p\ 2 \pi - p\ 2 p 2 - p 2 - ie K ' 

This is a standard integral that appears in the computation of radiative corrections to 



Coulomb scattering |T3[ and we solve it explicitly in the appendix. The leading non- 
analytic behavior of this integral is 

T 7T 6 |£| 5 



256 p) ' 

In order to Fourier transform this expression we need: 

A 



(3.5) 



/ 



d^e^lkl 5 = -2 15 3 3 5 2 ^. (3.6) 



r 14 ' 



Up to a constant the final result for the effective potential is 

(2) 225AgiV^ 

V eff<* 64 R5 M lS r U- V-<) 

This shows that this diagram indeed gives a contribution to the t> 6 /r 14 -term. To get 
the exact numerical coefficient one would have to evaluate all possible double graviton 



exchange diagrams. Moreover from the appendix we observe that the imaginary part of 
this amplitude is infrared divergent 0. This is familiar from the calculation of radiative 
corrections to Coulomb scattering in four-dimensions [|13]]. To cancel this divergence one 
gives a mass to the photon and the sum over all polarizations gives a finite amplitude. 
Polarizations for the probe graviton will have to be taken into account here as well at 
some point. This calculation seems much harder because graviton vertices are rather 
complicated. The easiest way to proceed might be to compute this four-graviton amplitude 
using type IIA string theory as we do in the next section. Of course, we expect the 
loop amplitude to reproduce the result of the classical calculation of ||. This is similar 
to the situation appearing in the four-dimensional Newton potential^ [|15]. Here general 
relativity corrections can be computed through graviton scattering amplitudes or by a 
classical calculation in which the zero component of the Schwarzschild metric is expanded. 

3.2. The Covariant Four- Graviton Amplitude 

It is curious to see that in the space-like compactification the general relativity correc- 
tion of the previous section can be computed from the covariant four-graviton amplitude 
0] |J . The sum of the tree level and one- loop contributions to the on-shell four-graviton 
amplitude in D dimensions is given by Q: 

kinematical factor x 2 T J - 1 t 2 T ^ V4 + ciO ( } + . . . , (3.8) 



where c\ is a calculable number determined by unitarity [[T6[ and s, t and u are the 
Mandelstam variables 



\2 



s = -(pi +P2) 

t = -( Pl +p 4 ) 2 , (3.9) 

U = -(pi +P3) 2 - 

Here 1 and 2 are incoming particles while 4 and 3 are outgoing and pi describe their corre- 
sponding momenta. The first term in formula ( |3.8|) represents the tree level contribution 



Infrared divergences in quantum gravity can be treated in the same manner as in quantum 
electrodynamics |14}| . 

We thank Steven Weinberg for pointing out this reference to us. 
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while the one-loop part of the amplitude g^ 1 ' is given by: 

9 {1) = K \I if^nr) [F 2 (a, r)] 10 "^ , (3.10) 

a' J (nary 

where a = ya! /R is a dimensionless parameter. The expression for F(r) can be found 
in Q . We are interested in the low-energy limit of this amplitude which corresponds to 
the zero slope limit a' — > 0. The one-loop contribution is divergent for D > 8 (here we 
consider D = 11). The finite part of (|3.8|) reduces to 

(kinematical factor) x I \- c\g ' + ...], (3-H) 

where 

g W = K 2 D c(-f) [J 7 (s, t) + 7 7 (t, s) + J 7 (a, u) + I 7 (u, s) + I y (t, u) + I y (u, t)] , (3.12) 

is the asymptotic value of g^ 1 ' and 7 = D/2 — 4. Here 0(7) is a constant determined by 
the dimensionality of the space-time 

1 /7T\7+l/2 , 

^ = 4(4) [sin(7r 7 )r(7+5/2)] \ (3.13) 

The function 1(7) is 



J 7 (M) = F +1 / \_ L t - (3-14) 



.! r 1 (i-xy +1 dx 

'0 sx — t(l — x) 
In .D = 11 we get 7 = 3/2 so that 0(7) is: 

^2 



c(3/2) = -^_, (3.15) 



and 6rJ 1} 



#0 = ~^4 K 11 [ J 3/2( S 5 *) + h/2& S) + 4/2(s, M) + h/2(u, s) + h/ 2 {t, u) + h/ 2 (u, tj\ . 

(3.16) 
The integrals ^3/2 can all be evaluated exactly but give rather complicated functions. 
Expanding around t = we are interested in the leading non-analyticQ term in t. It is 
given by: 

9P=i&l~+0{tV'). (3.17) 



The analytic terms correspond to contact terms after Fourier transforming. 



The kinematical factor in fl3.8| ) can be computed from the tree level contribution of this 
formula. We consider the case where the first graviton has vanishing transverse velocity, 
i.e. pi = (M, 0, M), where M = N s /R is the mass. The second incoming graviton has 
P2 = (E2,P2, to), where E 2 = p§ + w 2 is the relativistic energy, pi is the nine-dimensional 
transverse momentum and m = N p /R is the mass. To leading order in the velocity we 

have 

s = —Mmv 2 + . . . , 

t = -\k\ 2 + ..., (3.18) 

u = —s — t. 
Here we have set \j>2\ = rnv to leading order in v. As can be easily checked the result for 
the tree level contribution computed in |1|] and ||TI[ is reproduced with a kinematical factor 
— /■Ci 1 s 2 ii 2 . The tree level amplitude then becomes 

2 4 

A t , ee = K 2 11 S — = -K 2 11 M 2 m 2 ^-. (3.19) 

t \k\ 2 

In order to obtain the potential we are interested in the Fourier transformed of this ex- 
pression with respect to the momentum transfer. Using the relation 

(2tt) 9 |£|2 2(27r) 4 r 7 ' { ' 

we obtain 

1 f d 9 k ,^, 15 N n N^v 4 

In order to obtain this result we have included a normalization factor of l/\ / 2E for each 
state. 

Using the above kinematical factor we can obtain the form of the one-loop contribution 
to the potential. The one-loop amplitude is given by: 

Ai-ioop = -4iCi^iM 3 m 3 v 6 \k\ 5 . (3.22) 



v ^ = ^r I w'*'*- = is «?■ (3 - 21) 



We Fourier transform this expression using (|3.6| ). Up to a constant the result for the 
one-loop contribution to the potential is then: 



225 N±N : 



2 /V_ .,6 



Vi(r)ocJy p^t^T^^t- (3-23) 

y ' p 64 R 5 M}? r 14 v ; 



We observe that the covariant amplitude differs by a factor N p from the result calculated in 
using a classical picture. The difference is due to the fact that the covariant amplitude 
integrates over all eleven-dimensional momenta running in the loop while in the M(atrix) 
theory calculation only states with fixed p- are allowed to run through the loop. Both 
amplitudes are related^ by a phase factor times N p . This has to be taken into account 
when loop diagrams are calculated while tree diagrams naturally agree. 

4. Quantum Gravity Corrections 

M(atrix) theory not only allows us to compute corrections coming from general rela- 
tivity but the much harder quantum gravity corrections. They have, in general, fractional 
powers of k because a dimensionful parameter, the cutoff, is present in the calculation. 
An example of a correction of this sort is the v 8 /r 18 -term appearing in M(atrix) theory at 
two-loops. The result for this term in the effective potential can be easily obtained by a 
simple extension of the calculation performed in |TJJ 

_ 3 ■ 5 2 ■ 7 2 • 59 N?N P v 8 

2 11 i?7 M 24 r i8- ^ 

Susskind conjectured that this term is related to a four graviton amplitude involv- 
ing the i? 4 vertex of M-theory From the M-theory point of view the following is 
known. The leading term in the M-theory effective action is the conventional Einstein 
term. However there could be corrections coming from higher dimensional operators: 

R 2 + R s + R 4 + . . . . (4.2) 



Fradkin and Tseytlin [22| have shown that 

R m = for m = 0,1,2,3 (4.3) 

at one loop. The one-loop correction to the effective action takes then the form: 

S oc I dP-xy/^ghlyE?; (4.4) 



We thank David Gross and Lenny Susskind for discussions on this. 
Some discussion on this appeared in [18] [19] [20| [21]. 
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where An is a dimensionful cutoff. This cutoff can be determined from the type IIA 
string theory, as was done in |23| and [24|] . These authors decompactified the four-graviton 



amplitude of d = 9 closed superstring theory to calculate the R 4 interaction in eleven 
dimensions: 

S =^T I ^9 R + -573 f V^htsR 4 , (4.5) 

ZK 11 J Ki! J 



where c is a number that can be extracted from ||23|1 . It is plausible to use this result to 
compute the analog of (|4.1|) from the supergravity point of view. While it seems hard to do 
a precise calculation because we expect polarizations and fermions to be relevant, it is not 



too hard to see that the four graviton amplitude involving the vertex (|4.5|) indeed gives 



a contribution to ( fPD . We can use the probe-source picture of the previous sectional 
Higher-order local curvature invariants in eleven-dimensional supergravity do not change 



the form of /i__ pq| . The only effect of such a term is to change the action of the probe, 



so that the equation of motion for (p involves a new potential: 



2; \ 2 



d 2 _h 
dx 



n^r^ocf— =-) a T V + ... (4.6) 



This potential comes from the term (R- r - r ) 2 (R+i+j) 2 of the interaction vertex. The dots 
indicate other contributions in which we are not interested. Up to a constant the effective 
action from the supergravity calculation is: 

N 3 N 2 v 8 
S ^^WMft^ (4 ' 7) 

The dependences on», r N s , R and M p \ agree with (|4.7|), but we seem to find a 
disagreement for the dependence on N p . While the M(atrix) theory result is linear in N p , 
the M-theory result is proportional to Np. It seems plausible that the origin of this dis- 
crepancy is again the fact that in the M-theory calculation states with eleven-dimensional 
momentum are running through the loops, while in the M(atrix) theory calculation these 
states have a fixed p-. This changes the N p dependence of the result as we have seen in 
the previous sections. A more careful analysis including polarizations and fermions would 
have to be done to check if the numerical coefficient agrees. While the M(atrix) theory 



We thank Joe Polchinski for discussions on this. 
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calculation was an easy extension of our previous calculation [fT7| , the DLCQ M-theory 
counterpart seems to be much harder [26| . 



One may want to ask at this moment: why should the supergravity calculation agree 
with the M(atrix) theory calculation at all, especially because of the discrepancies found 
in the recent papers |27| |28| or why does the relativistic correction come out correctly in 
the calculation performed in []2|]? Probably the most natural explanations would be that 
the quantities computed in |27| and |28| receive corrections in the large iV limit, while the 
quantity computed in might obey a non-renormalization theorem that guarantees that 
the supergravity result is reproduced by M(atrix) theory for finite N. Our goal for the 
near future should then be to find out which quantities are not renormalized in the large 
N limit so that they can be computed with finite N M(atrix) theory. Maybe an argument 
along the lines of [E9j can be carried out here as well. 



Note Added 

These results have been presented in several seminars in October and November of this 
year and at the 3rd Workshop on Recent Developments in Theoretical Physics at CERN 
on December 9. 
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5. Appendix 

In this appendix we are going compute explicitly the integral: 

f 1 1 1 

/ - lim / dp y^ *l2\ l Z._,Z\2 ln 2_ p 2_ ie > ( 5A ) 



e ^° J \ Pf — P \ \ Pi — P \ P Pf 



We thank Lenny Susskind for discussions on this. 
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where p — \p\ and the same for pf. The above integrand can have poles when the denom- 
inator vanishes. In order to regularize them we introduce a "graviton mass" a. At the 
end of the calculation we will be interested in the leading order behavior in terms of a. 
Therefore we would like to compute the integral: 

111 



i" = lim / dp-— _ . ., ., . _, __, . ., ., ., 

e ^° J I Pf — P I + <T \Pi~ P \ + " P ~ Pf 

to leading order in a. First it is convenient to use the formula: 



ic 



(5.2) 



1 
ab 



1 



dz 



to write: 



where 



P~ Pf | 2 | P~ Pi I 2 



i (a(l + z)/2 + b(l-z)/2) 

i r 1 i d i 



2 ' 



(p-P) 2 + A 2 



dz, 



P=^[^ + z)p f + (l-z)p l 

A2 =a 2 + 1 (1 _ Z 2 )A ,2 5 



(5.3) 



(5.4) 



(5.5) 



and q = Pf ~ Pi is the momentum transfer. Using this identity I can be written in the 
form: 

• i 



1 
4,_, 



A<9A 



with 



/ = lim dp 



e^O 



1 



(P-P) 2 +A 2 P-PJ 



ie 



(5.6) 



We will compute the real and the imaginary part of the above integral separately. In order 
to do this we use the identity: 



1 



9 2 

P — Pf — It 



P.P. 



1 



9 2 

p z — pi 



+ in5(p 2 - p 2 f ) 



(5.7) 



where P.P. denotes the principal part and S is the Dirac delta function, we will first begin 
with the imaginary part. It is given by: 

8{p 2 -p)) 



Im/ = 7T / d p 



12 P 8 l v A2 + ( i 3_ ]?/) 2 



(5.8) 



Inserting in 



(p-P) 2 + A 2 
we can get the result for the imaginary part of I to leading order in a 



lml=0iog(^) + ^ 



(5.9) 
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The real part of / is 

Rel= / if'p— — ^-——;P-P- \ „■,. ' „■: ) • ( 5 - 10 ) 

to get 

Re/ = / dk I d 9 p " ^ (5.11) 

/o 7 (p-P) 2 + A 2 

This integral is straightforward to compute just a bit lengthy. Here we are only going to 

write down the terms of Re/ in which we are interested in. These are the terms that in 

the region for small momentum transfer are non-analytic. To leading order in a they are: 

Re/=— - t / f ox , |q| 5 I 1-^ + ^x 1 +... (5.12) 

4 {Ap 2 -q 2 f m I Q p 2 8Qp 4j V I 

The leading non-analytic behavior of / in terms of q is then given by: 

T n e \q\ 5 .n 5 p}, fa 2 \ ,_ 



256 pi 6 q 2 ' \q 
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